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4 Conclusions 



In the previous section we have presented results relevant for the weak field expansion of lat- 
tice gravity. We have shown the precise correspondence between continuum and lattice degrees of 
freedom within the context of such an expansion. To avoid unnecessary technical complications, 
most of our discussion has focused on the two-dimensional case, but the methods presented in 
this paper can be applied to any dimension. In the purely gravitational case, we have shown 
that the presence of zero modes is tied to the existence of a local invariance of the gravitational 
action. This invariance corresponds precisely to the diffeomorphisms in the continuum, with 
appropriate deformations in the edge lengths playing the role of local gauge variations of the 
metric in the continuum. 

We have then derived the Feynman rules for gravity coupled to a massless scalar field, to 
lowest order in the weak field expansion. Although the lattice Feynman rules for the edge length 
vertices appear to be rather unwieldy, they actually reduce to the familiar continuum form when 
re-expressed in terms of the weak field metric field, and we have presented this important result 
in detailed form. As an application, the two-dimensional conformal anomaly due to a massless 
scalar field was computed by diagrammatic methods. We have given explicitly the relevant 
lattice integrals, which involve among other things a new diagram, the tadpole term, which 
vanishes in the continuum but is necessary on the lattice for canceling unwanted terms. Finally 
we have shown that in the leading continuum approximation the expected continuum form for 
the anomaly is obtained, with the correct coefficient. 

The procedure followed here in deriving the Feynman rules for lattice gravity works in any 
dimension, including therefore the physical case of four dimensions. The lack of perturbative 
renormalizability in four dimensions is not ameliorated though by the presence of an explicit 
lattice cutoff, and non-perturbative searches for an ultraviolet fixed point and a lattice continuum 
limit are still required. 
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transformations (which change the local volumes and curvatures) in Eq. (3.20). The explicit 
form for the lattice diffeomorphisms, to lowest order in the lattice weak field expansion, is given 
in Eq. (3.18), which makes it obvious that such a decomposition can indeed be performed on 
the lattice. After rewriting the gravitational functional measure in terms of conformal and 
diffeomorphism degrees of freedom, 

[dg] = [dtp}[d X } [det(Z+Z)]5 , (3.98) 
which involves the Jacobian of the operator L, which in the continuum is determined from 

(L+L x )„ = V(V^ + V.x^ - Q^VXp) , (3-99) 

one has for the effective action contribution in the conformal gauge, and to lowest order in the 
weak field expansion, 

[det(i+i)]"2 ~ exp{-I eJJ (<p)} . (3.100) 

The lattice form of L depends on the specific form of the lattice gauge fixing term. On the 
lattice the functional integration is performed over the edge lengths, see Eq. (2.3). The lattice 
conformal gauge choice corresponds to an assignment of edge lengths such that 

*>■<»>= + + ® (3.1011 

The lattice fields <~p{n) are defined on the lattice vertices. So are the gauge degrees of freedom 
Xfi(n), as can be seen from Eq. (3.19). It should therefore be clear that the choice of lattice 
conformal gauge corresponds to a re-assignment of edge lengths about each vertex which leaves 
the local curvature unchanged but brings the induced metric into diagonal form; it corresponds 
to a choice of approximately right-angle triangles at each vertex [32]. 

A diagrammatic calculation, similar to the one for the scalar field contribution, gives in the 
continuum the celebrated result [15] 

13 1 

iWa(g) = (q»qv - <W 2 ) ^2 (^pftr - Spaq 2 ) • (3.102) 

As a consequence the total Liouville action for the field <p = -^R becomes 

Thus to lowest order in the weak field expansion the critical value of D for which the action 
vanishes is D c = 26, but this number is expected to get modified by higher order quantum 
corrections. In any case, for sufficiently large D one expects an instability to develop. Numer- 
ical nonperturbative studies of two-dimensional gravity suggest that in the lattice theory the 
correction is large, and one finds that the threshold of instability moves to D c ~ 13 [28]. It is 
unclear if this critical value can be regarded as truly universal, and independent for example on 
the detailed choice of gravitational measure. 
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which in real space becomes the well-known Liouville action 

Ieff(<p) = - / d 2 x [(d^ f + (A - A c ) e*] . (3.94) 

One therefore completely recovers the result derived perturbatively from the continuum Feynman 
rules, as given for example in [15]). 

In the continuum, the effective action term of Eq. (3.94) arising from the conformal anomaly 
can be written in covariant form as 

\ j d 2 x d 2 y R^(x) (x\—^—^\y) R^/g(y) , (3.95) 

where d 2 is the continuum covariant Laplacian, d 2 = d^y/gg^ d v \ and m 2 0. In two 
dimensions and for flat space, {x\-^j\y) ~ ^Mog|a; — y\. Using the correspondence between 
lattice and continuum curvature operators derived in [5], its lattice form gives rise to an effective 
long-range interaction between deficit angles of the type 

h Y S h \— o] h' , (3.96) 

hmgesft,/r 

where A is the nearest-neighbor covariant lattice Laplacian, as obtained from the discrete scalar 
action, and m 2 <— an infrared mass regulator. This result was given in [9]; see also the recent 
discussion in [35]. 

In general one cannot expect the lattice expression for the vacuum polarization to match 
completely the continuum expression already at the lowest order of perturbation theory. Since 
there is no small parameter controlling the weak field expansion in two dimensions, it is difficult 
to see why the first few orders on the lattice should suffice. Indeed in the weak field expansion 
the background lattice is usually taken to be regular, which leads to a set of somewhat preferred 
lattice directions for high momenta, which are close to the momentum cutoff at ±7r. It would 
seem though that this is an artifact of the choice of background lattice (which is necessarily rigid), 
and whose effects are eventually washed out when the fluctuations is the edge lengths are properly 
accounted for in higher order in the weak field expansion. The fact that the conformal mode 
in fact remains massless in two dimensions in the full numerical, non-perturbative treatment of 
the lattice theory was shown in [16,28], without the necessity of any sort of fine-tuning of bare 
parameters. 

The gravitational contribution to the effective action in the lattice conformal gauge can, 
at least in principle, be computed in a similar way. Let us sketch here how the analogous 
lattice calculation would proceed; a more detailed discussion of the relevant calculations will be 
presented elsewhere. In the continuum the metric perturbations are naturally decomposed into 
conformal and diffeomorphism parts, 

£>g^{x) = g^(x) Sip(x) + Vfj,Xu(x) + S7 u Xfj,( x ) ■ (3.97) 

where V„ denotes the covariant derivative. A similar decomposition can be done for the lattice 
degrees of freedom, by separating out the lattice gauge transformations [2] (which act on the 
vertices and change the edge lengths without changing the local curvatures) from the conformal 
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Analytical evaluations of the remaining integrals in terms of Bessel functions will be pre- 
sented elsewhere [34]. In the following we will calculate the graviton self-energy using the 
lattice Feynman rules in the leading continuum order, by doing a continuum approximation 
to the integrands valid for small momenta. The procedure is justified for integrands that are 
sharply peaked in the low momentum region, and neglects the effects due to presence of a high 
momentum cutoff. 

The calculation is most easily done by using the Feynman rules in terms of h^ v in the 
Lorentz covariant form, which were given before. The vacuum polarization loop contributing to 
the graviton self-energy then reduces to 

1 f d 2 p tnv(p,q)t a p(p,q) 



2 J (2vr) 2 p 2 (p+q^ 2 



t»v(p,q) = ^[6»vP-(p + q)-p»(pv + qv)-Pv(p» + q»)] ■ (3.86) 

The calculation of the integral is easily done using dimensional regularization. One obtains 

V^apiq) = {q 2 &nv - q»q») (? 2 <W - q a qp) ■ (3.87) 

As expected, the tadpole contribution T^^piq) to the graviton self-energy is 0. The remaining 
graviton self-energy contribution is given by 

^■iiv,ap(q) = V lll/}0 p(q) + T lll/}0 p(q) 

= (^ 2 <W - i^) ^ (9 2 <W - q a qp) • (3.88) 

The above result is valid for one real massless scalar field. For a _D-component scalar field, the 
above result would get multiplied by a factor of D. Then the effective action, obtained from 
integrating out the scalar degrees of freedom, and to lowest order in the weak field expansion, 
is given by 

i f d 2 q 

, W) 2 

with the scalar vacuum polarization given by 



In the lattice analog of the continuum conformal gauge, 

9^{x) = S^evW , (3.91) 
one can write for the scalar curvature 

{q^qu - ?>^q 2 )h^(q) = R(q) = q 2 <p(q) , (3.92) 



and therefore re-write the effective action in the form 

D [ d 2 q 

96W (2tt) 2 



lefffr) = 1 7S2 <p(q)qM-q) , (3-93) 
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T 



- r d 2 P qi r . 2 Pl . 

/ 7 TT: cos — < —4 sin — 

.J- v (2vr) 2 2 I V 2 • 

1 



4 sin — + 2 sm 



(3.81) 



The first three integrals are (^independent, while in the last three the g-dependence factorizes. 
One finds 



T, 



£2, £2 \ 



T, 



T 



T 



T, 



^2,^3 < 



3 J_ 

8 ~ 4^ 

£ _ _L 

8 ~ 4^ 

£ _ _L 

4 ~ 2^ 

ft + ft ^ f}_ 

2 7T 



COS 



ft! 



1 



cos(^) + - 

COs(f) + - 



(3.82) 



The terms in T tt)tj (q) are required to cancel some of the unwanted, including noncovariant , terms 
in V tt)tj (q). Also, a number of contributions to the lattice vacuum polarization can be shown to 
contribute to a renormalization of the cosmological constant. It should be noted here that the 
lattice form of the cosmological constant term (which just corresponds to the total area of the 
triangulated manifold) contains, in contrast to the continuum, momentum dependent terms [9]. 
The reason for this is that the lattice area terms couple neighboring edges, and lead therefore 
to some residual local interactions between the edges variables, as shown in Eq. (3.14). 

Eventually one needs to rotate the final answer for the vacuum polarization from the e 4 - to 
the variables, which is achieved to linear order (and for small q) by the matrix 




with e(q) 



V 

i i 

4 4 2 

V ■ h(q) (see Eq. (3.13)). Since the integration over the scalar gives 



(3.83) 



l eff 



d 2 q 



h(q) -11(h) -h(q) 



d 2 q 



e(q) ■ 11(e) • e(q) 



one obtains the correspondence between the two polarization quantities, 

11(h) = V' 1 ■ 11(e) -V , 
correct to lowest order in the weak field expansion. 



(3.84) 



(3.85) 
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(3.73) 



The g 2 -dependent terms are the ones that are of physical importance. Expanding the lattice 
integrals for small external momentum one obtains for example 

V tutl (ql = q 2 = q) = ±(l-^)+ ^-q 2 + 0(q A ) , (3.74) 

which can be compared with the continuum result for the metric vacuum polarization 

111111 = ^T f l 2 + ■ ( 3J5 ) 

907T 

A complete analytic evaluation of all the lattice integrals in the limit of small momentum is 
beyond the scope of this paper, and will be presented elsewhere. 

The expressions relevant for the tadpole diagram, written in component form, are 

T ei,M = I I \ 2 8111 ( — ) + 8111 ( — ) " o 8111 z 

(3.76) 



n (2vr) 2 I v 2 ' v 2 ' 2 v 2 

1 



4[E, sin 2 (^; 



^2,^2 (?) = / / TTTZyi sm (t) + 2 sm (t) - o sm 



x (2vr) 2 I v 2 7 v 2 7 2 v 2 

1 



4[E„ sm 2 (f) 



(3.77) 



x (2tt) 2 
1 



r e3je3 (g) = / / — — 2 ^ 3 sin 2 ( — ) + 3 sin 2 ( — ) - sin 



4[E„ «m 2 (^; 



(3.78) 



^ / ^ / i v d 2 p qi + q 2 f . 2/PU , o • 2/P2 



2/Pl +P2^] 1 



-2 sin 2 ( ^ 1 ' U x — — 5—5 — — . (3.79) 

1 2 V 4[£ At sin 2 (^)] 1 J 



1 



4[E, sin 2 (f; 



(3.80) 
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x r / \ i r r d 2 p r q2.p1. Pi + qi 

V (q) = — I I — < 2 cos — sm — sm 

3, 3W 2 J- v J- v (2 7r) 2 \ 2 2 2 

. qi . P2 . P2 + Q2 . pi + P2 . pi + qi + P2 + q2 ^ 2 

+2 cos — sm — sm 2 sm sm 

2 2 2 2 2 

1 (3.69) 



16 [Em ^ 2 (f)][E, sin 2 (^)] 



^i,e 2 (?) = o / / To — ^ (-^ sin — sm h cos — sm ■ 



£l ' £2W " 2 i_ T i_ T (2vr)n V 2 2 T ^ 2^ 2 

. Pi + ft + P2 + ft , , . P2 . P2 + ft . 

X sm X —2 sm — sm h 

2 ; v 2 2 

ft . Pi + P2 . Pi + ft + P2 + ft 

cos — sm sm 

2 2 2 

1 



16 [E, sin 2 (^)][E A1 sin 2 (^)]' 



x r / \ 1 r \( 1 • P 1 ■ P 1 + ^ , ft • P1+P2 

^ei^iftl = o / / 77^ — 7? i v sm — sm h cos — sm 



£l ' £3W - 2 i_ T i_ T (2vr)n V 2 2 T ^ 2^ 2 

s , . Pi + ft + P2 + ft s ^ , ft . Pi . Pi + ft . 

X sm X 2 cos — sm — sm h 

2 ; v 2 2 2 

ft . P2 . P2 + ft . Pi + P2 . Pi + ft + P2 + ft 

2 cos — sm — sm 2 sm sm 

2 2 2 2 2 

1 



16 [E„ sm^HE, sin 2 (^)]' 



(3.70) 



(3.71) 



t/ / \ i r r rf2 ^ // • P2 . P2 + ?2 . ft . P1+P2 

K E , e , (o) = - / / — < — 2 sin — sm h cos — sm 

2 ' 2 J-nJ-n (2vr) 2 I 2 2 2 2 

s , . Pi + ft + P2 + ft s ^ , q2 . pi . pi + ft . 

X sm X 2 cos — sm — sm h 

2 ; v 2 2 2 

ft . P2 . P2 + ft . Pi + P2 . Pi + ft + P2 + ft 

2 cos — sm — sm 2 sm sm 

2 2 2 2 2 

1 



" 16[E A1 sin 2 (^)][E A1 sin 2 (^)]- (3J2) 
For zero external momentum q all the quantities V^^.'s are infrared finite, and one finds 





= 


0) = 


1 

16 


0-1 




,e 2 (<l = 


0) = 


1 

16 


H 


Vc 3 


,e 3 (<l = 


0) = 


i( 
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3.2 Conformal Anomaly 



As an application, we will compute the graviton self-energy using the lattice Feynman rules 
developed above. There are two diagrams contributing to the lattice graviton self-energy, namely 
the vacuum polarization loop and the tadpole diagram shown in Figure 8. 





Fig. 8. Lowest order diagrams contributing to the conformal anomaly. 

The evaluation of the above diagrams corresponds to a functional integration over the scalar 
field, performed to lowest order in the weak field expansion, and with a scalar field measure 
deriving from the functional measure 



II dl^F £ [l]x 

edges ij 



It U d € ; 



(3.66) 



the specific form of the gravitational measure will not matter in the following calculation since 
only the scalar field is integrated over. The relevant expressions for the vacuum polarization 
loop diagram, written in component form, are then 



d 2 p 



• Pi ■ Pi 

-2 sin — sin 



qi 



(h ■ Pi + 

cos — sin 

2 2 



P2 



sin : 



2 

Pi + qi 



2 

P2 ■ 



Q2 



16 sin 2 (f)][E A1 sin 2 (^; 



(3.67) 



V fr 



1 r r d 2 p 

2 J- v J- v (2vr) 2 



. P2 ■ P2 + q2 

-2 sin — sin 

2 2 



qi . P1+P2 . pi + qi+ P2 + q2 

cos — sin sin 

2 2 2 

1 

16 sin 2 (f)][E A1 sin 2 (^)]' 



(3.68) 
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and its leading continuum approximation for small momenta is 

^(Pift: +P2ft) • (3.59) 
For the vertex h 22 (k)h 12 (k')cf)(p)(f)(q) ( = h 22 (k)h 21 (k')cf)(p)(f)(q) ), 

\ (- cos \ + 3 cos " 4 cos fcl+fc 2 2+fc ^ ) sin(f ) sin(f ) 

+ | (- cos ^ - 4 cos f + 3 cos sin (f ) sin(f ) 

+ cos(^) sin(a±^) sin(^±22) , (3.60) 
and its leading continuum approximation for small momenta is 

^(Pift: +P2ft) • (3.61) 

For the vertex h 12 (k)h 12 (k')4>(p)4>(q) ( = h 21 (k)h 21 (k r )4>(p)4>(q) = \ [h 12 (k)h 21 (k') 
(f>(p)<f>(<l) + h2i(k)h 12 (k')(f>( P )(f>(q)]), 

i(_ cos ^ + 3 cos^)sin(a) sin( a) 

+ \{- cos \ + 3 cos ^4^) sin (f ) sin(f ) , (3.62) 

and its leading continuum approximation for small momenta is 

\p-1 ■ (3-63) 

The above lattice Feynman rules for the 2-graviton 2-scalar vertex in terms of h^ v can be 
compared to the continuum Feynman rules, which are given by the following expression 

+ Ka {p^qp + ppq^) + {p^q a + p a q^) - ^ (p 

- 8 a p(Pv.qv + Pvq») + (SpvSap - ^iiaKp - S ll pS va )p ■ q] . (3.64) 

Again one sees that the 2-graviton 2-scalar vertex in the leading continuum approximation 
reduces completely to the continuum Feynman rules. 

In order to see how much the lattice Feynman rules depend on a particular lattice triangu- 
lation, one can derive the Feynman rules for the lattice obtained by a reflection with respect to 
the vertical axis (see Figure 6b). It turns out that these Feynman rules can simply be obtained 
by performing the following substitutions in all the above formulae for the vertices, 

Pi -pi,qi -qi,h -h,k[ -k[ , (3.65) 

which indeed corresponds to a reflection about the vertical axis. In this case the leading contin- 
uum approximation of the lattice Feynman rules are again identical to the continuum Feynman 
rules. One therefore concludes that the two inequivalent lattice triangulations give the same 
physical results, at least to the first leading continuum order (for momenta which are small 
compared to the ultraviolet cutoff). This should not come as a surprise, since the two lattices 
correspond to two equivalent parameterizations of flat space, with an action that is parameter- 
ization invariant, at least for small deformations. 
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for the vertex €2{k)e3(k l )<f)(p)<f)(q), 



_4 cos( fcl+fc 2 2+ ^ ) sin(^) sin(f ) - 4 cos(f) sin(f ) sin(f ) 

+ 2 cos(^) sin(£i±^) sin(ai±22 ). (3.51) 

After writing again the e 4 - fields in terms of the /i^j, fields (see Eq. (3.12)), one obtains the 
Feynman rules for the vertex hn(k)hn(k')(f>(p)(f>(q), 

(6 - \ cos \ + I cos _ 4 cos f ) sin(f ) sin(f ) 

+ (| cos ^ - 4 cos *i' + 2 cos fcl+fc2 + fcl ' +fc ^ - i cos f ) sin(f ) sin(f ) 

+ ( cos h. _ cos h+tk) sin(£i±^) sin(^±^) . (3.52) 

The above expression should be symmetrized in k and k' to reflect the fact that an edge is shared 
between two neighboring triangles. For small momenta its leading continuum approximation is 
given by 

- {-p • q + 4pigi) . (3.53) 

For the vertex h22{k)h22{k l )4>{p)4>{q), 

(3 cos h+tk _ 4 cos fcl+fc 2 2+fc ^ + 2 cos fcl+fc2 + fcl ' +fc ^ - \ cos \) sin(f) sm (f) 
+ (6 - 4 cos f + f cos " | cos t) sm (f) sin (f ) 

+ (cosf -cos sin(£i±^) sin(^±22) , (3.54) 

and its leading continuum approximation for small momenta is given by 

j (-P ■ q + 4p 2 ft:) • (3.55) 

For the vertex hn(k)h22(k')(f>(p)(f>(q), 

2(- cos f + cos + | cos ^ - cos *l±*2±*i - I cos f ) sin(f ) sin(f ) 

+ 2(- cos f + f cos ^ + cos ^ - cos fcl+fc 2 2+fc1 ' - i cos f ) sin(f ) sin(f ) 

+ (cos ^ - 2 cos + cos t) sm(a±^) sin(^±^) , (3.56) 

and its leading continuum approximation for small momenta is 

1 



p-q ■ (3.57) 



For the vertex hn(k)hi2{k')(t)(p)(t)(q) ( = h 11 (k)h 2 i(k')(f)(p)(f)(q) ), 

\ (- cos if - 4 cos f + 3 cos ^±^) sin(^) sin(f ) 
(- cos \ + 3 cos ^4^ " 4 cos fcl+fc 2 2+fcl ' ) sin(f ) sin(f ) 

+ cos(^) sin(£4^) sin(^±22) , (3.58) 



l 

2 



19 



and its leading continuum approximation is given by 



^ (-P2 42 +Pi qi) = ^ (p ■ q ~ 2p 2 ft:) • (3.42) 



For the vertex hi2{k)4>{p)4>{q) (= h2\{k)<f)(p)<f)(q) ) one has 



2 cos(f ) sin(f ) sin(f ) + 2 cos(Jf ) sin(f ) sin(f ) 

- 2 sin(^±^) sin(^±^) , (3.43) 

and its leading continuum approximation is given by 

-\{Piq2+P2qi) ■ (3.44) 

Written in Lorentz covariant form, the leading continuum approximations for the vertices 
h^u{k)(j)(p)(j)(q) can be grouped together as 

^(fi»vP- q- p»q» - p»q») (3.45) 

which is now identical to the usual continuum Feynman rule, derived from the original continuum 
action. 

One can proceed in a similar way for the higher order vertices. The Feynman rules for the 
2-scalar 2-graviton vertex written in the component forms are: 
for the vertex ei(k)ei(k')(f>(p)(f>(q), 

2 sin(f ) sin(f ) + cos( fcl+fc2 + fc i +k ? ) s in(f ) sin(f ) 

- \ cos(^±^) sin(£i±^) sin(*i±22 ) ; (3.46) 

for the vertex €2{k)€2{k')<f)(p)<f)(q), 

cos( fcl+fc2 + fcl ' +/ 4 ) + 2 

_ I cos (^±^) sin(Hi±^) sin(^) ; (3.47) 
for the vertex e^(k)e^(k')(j)(p)(j)(q), 

3 cos(^) sin(f ) sin(f ) + 3 cos(^) sin(f ) sin(f ) 

- sin(£i±^) sin(Si±^) ; (3.48) 



for the vertex €i(k)e2{k l )<f)(p)<f)(q), 



sin(f ) sin(f ) + 2 cos(^) sin(f ) sin(f ) 
- 2 cos(^^) sin(£i±^) sin(Si±^) ; (3.49) 



for the vertex ei(k)e^(k')(j)(jj)(j)(q) , 



k i) sin(f ) sin(f ) - 4 cos(^±^) sin(f ) sin(f ) 

+ 2 cos(^) sin(a±^) sin(^±22 ) ; (3.50) 
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Fig. 7. Labelling of momenta for the scalar- graviton vertices. 

The higher order terms give, after transforming to momentum space, the Feynman rules for 
the vertices. For the trilinear vertex associated with €i(k)<f)(p)<f)(q) one finds 

- 2 sm (y) sm (y) + cos (y) sm ( — ^ — ) sm ( — ^ — ) ; ( 3 - 34 ) 
for the vertex e 2 (k)<f)(p)<f)(q) one obtains 

— 2 sin — sin — + cos — sin sin ; 3.35 



2' v 2 ' v 2 
and for the vertex e^(k)(j)(p)(j)(q) one has 



h.) sin( a) sin(f ) + 2 cos( ^) sin(f ) sin(f ) 

- 2 sinfSii^) sin(Si±^) . (3.36) 



In order to compare the lattice Feynman rules with the usual continuum ones, one needs to 
perform a transformation from the e 4 - variables to the metric fluctuation h^ v . The correspondence 
between the two is given in Eqs. (2.13), (3.12) and (3.13). However, one must be careful in 
doing the transformation since 63 contains contributions to all fin, h 22 , an( i ^12- The exact 
correspondence is given by the following relation. After writing 

ae 1 + be 2 + ce 3 + Ae\ + Be\ + Ce\ 

+ De x e 2 + £ ei e 3 + Fe 2 e 3 , (3.37) 

one can use the relationships between e 4 - and h^ v of Eq. (3.12) to obtain 



(f + |) h u + (| + I ) h 22 + \ h 12 + \ h 



21 



A , C_ 1 E\ h 2 1 ( b_ _c_ _j_ B_ _j_ C_ _j_ F_\ ^2 

"8 32 T 4 T 16 T 8 i ""H V 8 32 i "4 i "l6 i "8 / ""22 

16 1 ¥ + T + f + f ) hl1 ^ 22 + (~T& + ¥ + f ) ^11 ^12 
C , F\ u U , f c , C , f> 



'll + ¥ + 8") ^ n ^21 + (~JQ + 8" + 8" ) ^22 h 12 

-fe + § + f ) h 22 h 21 + (-£ + H) ^ + (__£. + |) /, 12 /, 21 . (3.38) 



With the aid of the above equation, one can then easily rewrite the Feynman rules in terms of 
h^ v (\? . For the vertex hn(k)<f>(p)<f>(q) one has 

-4(1-1 cos h_ } sm( a } sm( ^ } + 2 cos( ^ } sin( ^ } sin(f } 

-2(l-cos^) sin(£i±^) sin(^±^) . (3.39) 
When expanded out for small momenta it gives 

\{-Pi1i +P2ft:) = ^(p-q-2p 1 q 1 ) . (3.40) 
For the vertex h 22 (k)<f)(p)<f)(q) one obtains 

2 cos(4f ) sin(^) sin(f ) - 4 (1 - \ cos \) sin(f ) sin(f ) 

-2(l-cos^) sin(£i|^) sin(^±^) , (3.41) 
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following expressions 



^ 1 1 



7 — 7T e l ~r 7T e 3 



(1 + ei) 2 4 2 1 ' 2 



1 o 1 o 3 Q 

2 1 4 2 4 3 

le 1 e 2 -e 1 e 3 -e 2 e 3 + 0(ef) ; (3.30) 



4 } 1_1 ,1 

(l + e 2 ) 2 ' 4 2 62 + 2 63 



1 9 1 9 3 9 

- ef + - e 2 , + - d 

4 1 2 2 4 3 

^e 1 e 2 -e 1 e 3 -e 2 e 3 + 0(ef) ; (3.31) 



4? 1 1 1 



T e l + 7 e 2 - 77 e 3 



(l + e 3 ) 2 4 1 4 z 2 

- - e 2 - - d - - d 

8 1 8 2 4 3 



2 2 2 

The 0-th order term in e gives the scalar field propagator 

1 



1 e 1 e 2 + le 1 e 3 + le 2 e 3 + 0(ef) . (3.32) 



4E>n 2 (W 



(3.33) 



which is the usual scalar propagator for the square lattice, and coincides with the continuum 
one for small momenta. 
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(3.27) 



while we still define the scalar field on the vertices, and therefore the Fourier transform in the 
usual way, 

Hn) = f f e~ tP ' n HP) ■ (3-28) 



■J-» (2vr) 2 

These formulae are easy to generalize to higher dimensions when a simplicial subdivision of a 
hypercubic lattice is employed, as first suggested in [3]. 




Fig. 6. Fwo equivalent triangulation of flat space, based on different subdivisions of the square 
lattice. 

The kinetic energy term for the scalar field can naturally be decomposed as a sum of three 
terms 

<ij> 13 



, ,, 13 1 + ei 7 2 ^ ,J M + e 2 

horizontal vertical 



2 

, 1 AWffo-fov 



\ E ^(y^)' - (3.29) 

diagonal 



where < ij > labels an edge connecting sites i and j. We have separated the sum over all edges 
into sums over horizontal, vertical, and diagonal edges. The series expansion of each term in 
the sums with respect to an edge length fluctuation in a particular direction is given by the 
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where the x„'s represent continuous gauge transformations defined on the lattice vertices [33]. 
The invariance of the quantum theory is again broken by the triangle inequalities, which in 
one dimension reduces to the requirement that the edge lengths be positive. Such a breaking 
is unavoidable in any lattice regularization as it cannot preserve the invariance under scale 
transformations, which are just special cases of diffeomorphisms. 

The weak field expansion for the purely gravitational part can be carried out to higher order, 
and the Feynman rules for the vertices of order h 3 , h 4 , ... in the i? 2 -action of Eq. (2.12) can be 
derived. Since their expressions are rather complicated, they will not be recorded here. 

3.1 Feynman Rules 

Let us consider next the scalar action of Eq. (2.21), 

\ I d 2 x^-gg^d^dA ~ \ A tJ (^M 2 . (3.23) 



In the continuum, the Feynman rules are obtained by expanding out the action in the weak 
fields h^v^x), 

\ j d 2 x^r d^d v <t> = \ j d 2 x(d<t>) 2 + \ j d 2 x h^{\s^{d4>? - d^d v <t>} + o(h 2 ) , 

(3.24) 

and by then transforming the resulting expressions to momentum space. 

On the lattice the action is again expanded in the small fluctuation fields e 4 -, which depend 
on the specific choice of parameterization for the flat background lattice. We will first discuss 
the lattice of Figure 6a, with the choice of labelling according to Figure 5. Figures 6a and 
6b correspond to two different gauge choices for the background metric which are physically 
equivalent; there are many others. The fluctuations in the edge lengths e 4 - (see Eq. (3.3)) and 
the scalar fields (f> at the point i, j steps in one coordinate direction and k steps in the other 
coordinate direction from the origin, are related to the corresponding e 4 - and (f> at the origin by 

e? + *> =-M4 0) , (3-25) 

where u>i = e~ tk ' and k{ is the momentum in the direction i. 

fU+k) = JWUf) , (3-26) 

where uj'i = e~ vp * . In practice it is actually more convenient to redefine the edge variables at 
the midpoints of the links, since this choice neatly removes later a set of complex phase factors 
from the Feynman rules. For the edge lengths we therefore define the lattice Fourier transforms 
as 
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for operator averages, and such a term is not needed, as in ordinary lattice formulations of gauge 
theories. 

The above zero modes correspond to the lattice analogs of diffeomorphisms. It is easy to see 
that the eigenvectors corresponding to the two zero modes can be written as 




/ 1 - wi 

l-u 2 



X2M 



(3.18) 



where Xi( w ) an( i X2( u ) are arbitrary functions (the above result is not restricted to two dimen- 
sions; completely analogous zero modes are found for the Regge action in three [13] and four 
[3] dimensions, leading to expressions rather similar to Eq. (3.18), with d gauge zero modes in 
d dimensions; their explicit form can be found in the quoted references). In position space one 
then has 



ei(ra) = Xi(n) - Xi(n + £1) 
£2{n) = X2{n) - X2(n + fi 2 ) 

^(n) = \xi{n) + \x2{n) - \xi{n + £1 + £2) - \xi{n + £1 + £2) 

(3.19) 

Using the correct relation between induced metric perturbations and edge length variations, 

_ ( sq \(6ii - sq - sq) 
911 ~ { Wl - sq - sq } sq 

one can easily show that the above corresponds to the discrete analog of the familiar expression 

Sg^ = - d^Xv - duXfi , (3.21) 

which indeed describes the correct gauge variation in the weak field limit. In the discrete 
case it reflects the invariance of the lattice action under local deformations of the simplicial 
manifold which leave the local curvatures unchanged [2]. The above relationships express in the 
continuum the well-known fact that metrics related by a coordinate transformation describe the 
same physical manifold. Since the continuum metric degrees of freedom correspond on the lattice 
to the values of edge lengths squared, one expects to find analogous deformations of the edge 
lengths that leave the lattice geometry invariant, the latter being specified by the local lattice 
areas and curvatures, in accordance with the principle of discussing the geometric properties of 
the lattice theory in terms of lattice quantities only. This invariance is spoiled by the presence 
of the triangle inequalities, which places a constraint on how far the individual edge lengths can 
be deformed. In the perturbative, weak field expansion about a fixed background the triangle 
inequalities are not seen to any order in perturbation theory, they represent non-perturbative 
constraints. 

The above invariance of the lattice action is a less trivial realization of the exact local gauge 
invariance found already in one dimension 

Sl n = Xn+l ~ Xn , (3.22) 



(3.20) 
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which can be inverted to give 



e 2 = ^h 2 2 - \h\ 2 + 0(h 22 ) 

e 3 = + h 22 + 2h 12 ) - i(/i u + h 22 + 2h 12 f + 0(h 3 ) 

(3.13) 

which was used in Eq. (3.7). Thus the matrix M w was brought in the continuum form after 
performing a suitable local rotation from the local edge lengths to the local metric components. 

A similar weak field expansion can be performed for the cosmological constant term, although 
strictly speaking flat space is no longer a classical solution in the presence of such a term [9]. 
One then obtains a contribution to the second variation of the action of the form 



(3.14) 



In the weak field limit, and with the same change of variables as described for the matrix M w , 
this leads to 

/ 1 -1 \ 

-110 +0(k) . (3.15) 
y 4 / 

The local gauge invariance of the i? 2 -action is reflected in the presence of two exact zero 
modes in M w of Eq. (3.6). As discussed in Ref. [32], the eigenvalues of the matrix M w are 
given by 

Ax = 
A 2 = 

A 3 = 24 - 9(^i + uji + tjj 2 + uj 2 ) + 4(^1^2 + ^1^2) 




Li = - x - 

2 



-U\u\ + Lo\l0 2 + LO\ Lo\ + Lo\l0 2 



(3.16) 



It should be noted that the exact zero modes appear for arbitrary Ui, and not just for small 
momenta. 

It is clear that if one were interested in doing lattice perturbation theory with such an R 2 
action, one would have to add a lattice gauge fixing term to remove the zero modes, such as the 
lattice analog of 

^fe 1 ^[i)<r) 2 , (3.17) 

and then add the necessary Fadeev-Popov nonlocal ghost determinant. A similar term would 
have to be included as well if one were to pick the lattice analog of the conformal gauge [14]. 
If one is not doing perturbation theory, then of course the contribution of the zero modes will 
cancel out between the numerator and denominator in the Feynman path integral representation 



12 



The properties of M 8 j are best studied by going to momentum space. One assumes that the 
fluctuation e 4 - at the point i, j steps in one coordinate direction and k steps in the other coordinate 
direction from the origin, is related to the corresponding e 4 - at the origin by 

e? + *> = ^ 2 M 0) , (3.5) 

where u>i = e~ tk ' and k{ is the momentum in the direction i. The matrix M then reduces to a 
3x3 matrix M w with components given by [9] 

(AL)n = 2 + - 2u> 2 - 2wiw 2 + ^1^2 + c - c - 

(Muj)i2 = 2 - - UJ 2 - ^1^2 ~ ^1^2 ~ Ojj - Ljj + ljj\uJ 2 + + ^1^2 

(Muj)i3 = 2(-l + 2^1 - U! 1 + U) 2 ~ ^2 ~ ^1^2 + 2bJiUJ 2 +U 2 - UJlbjl ~ UiUJ 2 ) 

{M^)^ = 4(2 - 2^i - 2uj 2 + ^1^2 + ^1^2 + c.c.) 



(3.6) 



with the other components easily obtained by symmetry. The change of variables 



e' x = ei e' 2 = e 2 e' 3 = ±(ei + e 2 ) + e 3 • (3.7) 
leads for small momenta to the matrix M' w given by 

M L = /4 k i k 2 k t -2k\k 2 +0(k 5 ) . (3.8) 
\ -2km -2/^2 4kj k 2 ) 

This expression is identical to what one obtains from the corresponding weak-field limit in the 
continuum theory. To see this, one defines as usual the small fluctuation field h^ v about flat 
space, which then gives 

VgR 2 = (V22 + V11 - 2V12) 2 + 0(/* 3 ) • (3-9) 

In momentum space, each derivative d v produces a factor of k u , and one has 

^JgR 2 = y Fi(W hpa- , (3.10) 

where V^pa- coincides with M' above (when the metric components are re-labelled according to 
11 1, 22 2, 12 3). 

It is easy to see the reason for the change of variables in Eq. (3.7). Given the three edges 
in Figure 5, one can write for the metric at the origin 

^-\\{ti-i\-il) q ) • (3 - n) 

Inserting l{ = /°(1 + e 8 ), with /° = 1 for the body principals (i = 1,2) and /° = y/2 for the 
diagonal (i = 3), one then obtains 

^2 



/^(1 + ex) 2 = 1 + hu 
\\ = (1 + e 2 ) 2 = 1 + h 22 
ll 2 3 = (l + e 3 ) 2 = 1 + ±(/in + h 22 ) + h u 



(3.12) 



11 



As a result, the existence of gravitational waves and gravitons in the discrete lattice theory has 
been established (indeed it is the only lattice theory of gravity for which such a result has been 
obtained). 

In the following we shall consider in detail only the two-dimensional case, although similar 
calculations can in principle be performed in higher dimensions, with considerable more algebraic 
effort. In pure gravity case the Einstein- Regge action is a topological invariant in two dimensions, 
and one has to consider the next non-trivial invariant contribution to the action. We shall 
therefore consider a two-dimensional lattice with the higher derivative action of Eq. (2.12) and 
A = 0, 

1(f ) = 4a ^ ■ (3-2) 

hinges h ^ 

The weak field expansion for such a term has largely been done in [9], and we will first recall 
here the main results. Since flat space is a classical solution for such an R 2 — type action, one 
can take as a background space a network of unit squares divided into triangles by drawing in 
parallel sets of diagonals (see Figure 5). This is one of an infinite number of possible choices 
for the background lattice, and a rather convenient one. Physical results should in the end be 
insensitive to the choice of the background lattice used as a starting point for the weak field 
expansion. 




Fig. 5. Notation for the weak-field expansion about the rigid square lattice. 

It is also convenient to use the binary notation for vertices described in references [3]. As 
discussed in the previous section, the edge lengths on the lattice correspond to the metric degrees 
of freedom in the continuum. The edge lengths are thus allowed to fluctuate around their flat 
space values, 

= + e,-) , (3.3) 

with /° = I® = 1 and /§ = y/2 for our choice of lattice. The second variation of the action is then 
expressed as a quadratic form in the e's, 

S 2 I = 4aJ2^M tJ e J . (3.4) 

ij 



10 



with hi is the length of the edge dual to l\. The baricentric dihedral volume for the same edge 
would be simply 

A h = (A 123 + A 134 )/3 • (2.23) 

It is well known that one of the disadvantages of the Voronoi construction is the lack of positivity 
of the dual volumes, as already pointed out in [5]. Thus some of the weights appearing in 
Eq. (2.21) can be negative for such an action. On the other hand, for the baricentric subdivision 
this problem does not arise, as the areas Aij are always positive due to the enforcement of the 
triangle inequalities. It is immediate to generalize the action of Eq. (2.21) to higher dimensions, 
with the two-dimensional Voronoi volumes replaced by their higher dimensional analogs. 

Mass and curvature terms can be added to the action, so that the total scalar action contri- 
bution becomes 

I(l\4 >) = \Y,^{^^) 2 +\Y,^(ra 2 + iRi)^ . (2.24) 

<ij> lJ i 

The term containing the discrete analog of the scalar curvature involves the quantity 

AiRi = Y,6h~VgR ■ (2.25) 

hDi 

In the above expression for the scalar action, Aij is the area associated with the edge Z 4 -j, while 
Ai is associated with the site i. Again there is more than one way to define the volume element 
Ai, [5], but under reasonable assumptions, such as positivity, one expects to get equivalent 
results in the lattice continuum limit. In the following we shall only consider the simplest form 
for the scalar action, with m 2 = £ = 0. 

3 Lattice Weak field Expansion 

One of the simplest problems which can be studied analytically in the continuum as well as 
on the lattice is the analysis of small fluctuations about some classical background solution. In 
the continuum, the weak field expansion is often performed by expanding the metric and the 
action about flat Euclidean space 

9iiu = <W + k/ V • ( 3 - 1 ) 

In four dimensions k = ^/32ttG, which shows that the weak field expansion there corresponds 
to an expansion in powers of G. In two dimensions this is no longer the case and the relation 
between k and G is lost; instead one should regard k as a dimensionless expansion parameter 
which is eventually set to one, k = 1, at the end of the calculation. The procedure will be 
sensible as long as wildly fluctuating geometries are not important in two dimensions (on the 
lattice or in the continuum). The influence of the latter configurations can only be studied by 
numerical simulations of the full path integral [9,16]. 

In the lattice case the weak field calculations can be carried out in three [13] and four [3] 
dimensional flat background space with the Regge- Einstein action. One finds that the Regge 
gravity propagator indeed agrees exactly with the continuum result [31] in the weak-field limit. 
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For the scalar field derivatives one writes [29,30] 



which corresponds to introducing finite lattice differences defined in the usual way by 

8^ (A^) 8 = <f> i+ll - <f>i . (2.18) 

Here the index fj, labels the possible directions in which one can move from a point in a given 
triangle. The discrete scalar field action then takes the form 

J ( /2 ^) = Ye E T- K 2 (^ - <h)(<h - <h) + ilifo - <h){<h - <h) + il(<h - - <hj\ ■ 

A 

(2-19) 

Using the identity 

(fr - <t>j){<t>i -<t>k) = \ \((f>i ~ <t>j? + (fr ~ ~ ((f>j ~ <t>kf\ , (2-20) 



^^iE^fM 1 ) 2 , (2-21) 



one obtains after some re- arrangements the simpler expression [29] 

<ij> ■> 

where Aij is the dual (Voronoi) area associated with the edge ij. In terms of the edge length 
lij and the dual edge length hij, connecting neighboring vertices in the dual lattice, one has 
Aij = \hijUj (see Figure 4). Other choices for the lattice subdivision will lead to a similar 
formula for the lattice action, but with the Voronoi dual volumes replaced by their appropriate 
counterparts in the new lattice. 




Fig. 4- Dual area associated with the edge l\ (shaded area), and the corresponding dual link h\. 
For the edge of length 1 the dihedral dual volume contribution is given by 

16A 123 16A 134 2 V 
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2.2 Scalar Field 



A scalar field can be introduced, as the simplest type of dynamical matter that can be 
coupled to the gravitational degrees of freedom. The continuum action is 

I\g, ( f>]=±Jd 2 x y /g[g^d»(f>d v (f>+(m 2 + ZR)(f> 2 ] , (2.13) 

The dimensionless coupling £ is arbitrary. Two special cases are the minimal (£ = 0) and the 
conformal (£ = |) coupling case; in the following we will mostly consider the case £ = 0. 

On the lattice consider a scalar 4>i and define this field at the vertices of the simplices. The 
corresponding lattice action can be obtained through the usual procedure which replaces the 
original continuum metric with the induced metric on the lattice, written in terms of the edge 
lengths [3,13]. Here we shall consider only the two-dimensional case; the generalization to higher 
dimensions is straightforward. It is convenient to use the notation of Figure 3, which will bring 
out more readily the symmetries of the resulting lattice action. Here coordinates will be picked 
in each triangle along the (1,2) and (1,3) directions. 




Fig. 3. Labelling of edges and scalar fields used in the construction of the scalar field action. 

To construct the scalar lattice action, one performs in two dimensions the replacement 

(l<-/? + /3 + <3) ¥ -' } f + q> ) ■ P.") 

which then gives 

det g^(x) det 9ij (A) = \ {2(/ 1 2 / 2 2 + l\l\ + iffi) -l\-l\- l 4 3 } = AA\ , (2.15) 
and also 
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and again join the resulting vertices. The baricentric dihedral volume is given simply by 

A d (l 2 ) = A/3 . (2.8) 
For the baricentric subdivision one then has simply 

Ah = l £ A t 

triangles t 
meeting at h 

Ah can also be taken to be the area of the cell surrounding h in the dual lattice, with 

Ah = J2 A 'i > ( 2 - 10 ) 

triangles t 
meeting at h 

with the dual area contribution for each triangle A,i given in Eq. (2.7). In general, if the original 
lattice has local coordination number qi at the site i, then the dual cell centered on i will have 
qi faces. A fairly complete set of formulae for dual volumes relevant for lattice gravity and their 
derivation can be found in [5]. In the following we shall refer to the Voronoi cell construction 
as the "dual subdivision", while we will call the baricentric cell construction the "baricentric 
subdivision". 

It is well known that two-dimensional Einstein gravity is trivial because the Einstein action 
is constant and the Ricci tensor vanishes identically. When a cosmological constant term and a 
curvature-squared term are included in the action, 

1 = J d 2 x^^X-kR + aR 2 ^ , (2.11) 

the classical solutions have constant curvature with R = ±^/A/a (there being no real solutions 
for A < 0). The theory with the Einstein action and a cosmological constant is metrically 
trivial, having neither dynamical degrees of freedom nor field equations, although non-trivial 
interactions can arise from the functional measure. The lattice action corresponding to pure 
gravity is 

I(l 2 ) = Y J A h [\-2kR h + aR 2 h ], (2.12) 

h 

with local volume element Ah and the local curvature given by Rh = 2Sh/Ah- In the limit 
of small fluctuations around a smooth background, I{1 2 ) corresponds to the above continuum 
action [9]. For a manifold of fixed topology the term proportional to k can be dropped, since 
J2h fih = 27rx, where x 1S the Euler characteristic. The curvature-squared leads to non-trivial 
interactions in two dimensions, although the resulting theory is not unitary. 

A number of results have been obtained from the above pure gravitational action. Argu- 
ments based on perturbation theory about two dimensions (where the gravitational coupling 
is dimensionless and the Einstein theory becomes renormalizable) suggest that there should be 
no non-trivial ultraviolet fixed point of the renormalization group in two dimensions. Explicit 
calculations in the lattice theory have shown conclusively that this is indeed the case in the 
absence of matter [9,16,10], as well as in the presence of scalar matter for a sufficiently small 
number of components [28]. 



(2.9) 
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Fig. 1. Dual area A,i associated with vertex 0, and the corresponding dihedral angle 

It is useful to introduce a dual lattice following, for example, the Dirichlet-Voronoi cell 
construction (see [26,27] and references therein), which consists in introducing perpendicular 
bisectors in each triangle and joining the resulting vertices. It provides for a natural subdivision 
of the original lattice in a set of non-overlapping exhaustive cells (see Figure 2.), and furthermore 
has a natural generalization to higher dimensions. 




Fig. 2. Original simplicial lattice (continuous lines) and dual lattice (dotted lines) in two dimensions. 
The shaded region corresponds to the dual area associated with vertex 0. 

The vertices of the original lattice then reside on circumscribed circles, centered on the 
vertices of the dual lattice. For the vertex the dihedral dual volume contribution is given by 




It should be pointed out that the above subdivision is not unique. Alternatively, one can 
introduce a baricenter for each triangle, defined as the point equidistant from all three vertices, 
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zero modes [3] (the triangle inequalities, which violate gauge invariance by imposing a cutoff 

on the gauge orbits, are in fact not seen to any order in the weak field expansion). + A gauge 

fixing term is not required in non-perturbative studies, since the contributions from the gauge 

zero modes are expected to cancel exactly between numerator and denominator in the functional 

integral for observables, as discussed for example in [4]. 

For a given action, the dynamics of the lattice will give rise to some average lattice spacing 
1 

a = [(/ 2 )] 2 , which in turn will supply the ultraviolet cutoff needed to define the quantum 
theory. It should be stressed that in the following we shall restrict our attention to the lattice 
theory, which is defined in terms of its lattice degrees of freedom only. Since it is our purpose to 
describe an ultraviolet regulated theory of quantum gravity, we shall follow the usual procedure 
adopted when discussing lattice field theories, and describe, in the spirit of Regge's original 
idea, the model exclusively in terms of its primary, lattice degrees of freedom: the squared 
edge lengths. As such, the theory does not require any additional ad-hoc regulators for defining 
conical singularities, for example. 

2.1 Curvature and Discretized Action 



In simplicial gravity the curvature is concentrated on the hinges, which are subspaces of 
dimensions d — 2, and is entirely determined from the assignment of the edge lengths. In two 
dimensions the hinges correspond to the vertices, and 8h, the deficit angle at a hinge, is defined 

by 

6 h = 2ir- °t ' ( 2 - 4 ) 



triangles t 
meeting at h 



where 9 t is the dihedral angle associated with the triangle / at the vertex h (see Figure 1). In d 
dimensions several rf-simplices meet on a (d — 2)-dimensional hinge, with the deficit angle defined 

by 

6 h (l 2 ) = 2vr - J2 ^( /2 ) • ( 2 - 5 ) 



d — simplices 
meeting on h 



Useful formulas for the cosine of the dihedral angles can be found in [5]. In two dimensions the 
dihedral angle is obtained from 



72 , 72 _ 72 

cos^ = /oi +; o2 / /l2 . (2.6) 

2«01«02 



(for the labelling see Figure 1). 



+ A small breaking of gauge invariance causes well known problems in perturbation theory. There are 
rather convincing arguments that this is not necessarily the case non-perturbatively, if the breaking can 
be considered small [24,25], as in the present case 
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correspondence between the lattice and continuum theories, and bring out the role of local gauge 
invariance in the lattice theory. We shall then discuss the Feynman rules for gravity coupled to 
a scalar field, and as an application compute the conformal anomaly in two dimensions. The 
final section contains our conclusions. 



2 The Discretized Theory 



In this section we shall briefly review the construction of the action describing the gravita- 
tional field on the lattice, and use the occasion to define the notation used later in the paper. 
In concrete examples we will often refer, because of its simplicity, to the two-dimensional case, 
where a number of results can be derived easily and transparently. In a number of instances 
though important aspects of the discussion will be quite general, and not restricted to specific 
aspects of the two-dimensional case. 

It is well known that in two dimensions quantum gravity can be defined on a two-dimensional 
surface consisting of a network of flat triangles. The underlying lattice may be constructed in 
a number of ways. Points may be distributed randomly on the surface and then joined to 
form triangles according to some algorithm. An alternative procedure is to start with a regular 
lattice, like a tessellation of the two sphere or a lattice of squares divided into triangles by 
drawing in parallel sets of diagonals, and then allow the edge lengths to vary, which will introduce 
curvature localized on the vertices. For arbitrary assignments of edge lengths, consistent with the 
imposition of the triangle inequalities constraints, such a lattice is in general far from regular, and 
resembles more a random lattice. In the following though we will narrow down the discussion, 
and think of the "regular" lattice as consisting of a network of triangles with a fixed coordination 
number of six, although many of the results in this work are expected to be quite general and 
should not depend significantly on the specific choice of local coordination numbers. 

The elementary degrees of freedom on the lattice are the edge lengths, with the correspon- 
dence between continuum and lattice degrees of freedom given locally by 

{9 xcM - {^} l=1 ... Nl > 

where the index i ranges over all N\ edges in the lattice. From the well known relationship 
between the induced metric in a simplex and its squared edge lengths, 

9iiV 2 ) = + ■ (2-2) 

one then has the essentially unique functional measure contribution 

/*oo 

Y[dg^{x) -+ / Y[dl* , (2.3) 

x J ° i 

supplemented by the additional constraint that the triangle inequalities be satisfied for all quan- 
tum fluctuations of the edge lengths [2]. Both the continuum metric and the lattice edge 
lengths imply some redundancy due to local gauge invariance of the action, which therefore re- 
quires gauge-fixing when performing perturbation theory due to the presence of the exact gauge 
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1 Introduction 



In the quantization of gravitational interactions one expects non-perturbative effects to play 
an important role. One formulation available for studying such effects is Regge's simplicial 
lattice theory of gravity [I]. It is the only lattice model with a local gauge invariance [2], and 
the only model known to contain gravitons in four dimensions [3]. A number of fundamental 
issues in quantum gravity, such as the existence of non-trivial ultraviolet fixed points of the 
renormalization group in four dimensions and the recovery of general relativity at large distances, 
can in principle be addressed in such a model. The presence of a local gauge invariance, which 
is analogous to the diffeomorphism group in the continuum, makes the model attractive as a 
regulated theory of gravity [4], while the existence of a phase transition in three and four 
dimensions [5,6,7,8] (but not in two [9,10]) suggests the existence of a (somewhat unusual) 
lattice continuum limit. A detailed discussion of the properties of the two phases characterizing 
four- dimensional gravity, and of the associated critical exponents, can be found in [7]. Recently 
calculations have progressed to the point that a first calculation of the Newtonian potential from 
the correlation of heavy particle world lines, following the proposal of [11], has become feasible 
[12]. These results indicate that in the lattice quantum theory the potential is indeed attractive, 
and has the correct heavy mass dependence. In the same work a general scaling theory for 
gravitational correlations, valid in the vicinity of the fixed point, was put forward. 

In view of this recent progress it would appear desirable to further elucidate the corre- 
spondence between continuum and lattice theories. The weak field expansion is available to 
systematically develop this correspondence, and it is well known that such an expansion can be 
carried out in both formulations. Not unexpectedly, it is technically somewhat more complex in 
the lattice theory due to the presence of additional vertices, as happens in ordinary lattice gauge 
theories. In the past most perturbative studies of lattice gravity have focused on the lowest order 
terms, and in particular the lattice graviton propagators [3,9,13]. As such, these did not probe 
directly important, genuinely quantum-mechanical, aspects of the theory. A systematic weak 
field expansion is generally useful, since it allows one in principle to determine subleading lattice 
corrections to the continuum results, which can be relevant in the analysis of the numerical 
non-perturbative results in the full theory. 

More importantly, the weak field expansion can be used to compare with known results in 
the continuum, and some are known in two dimensions [14,15]. A related motivation comes 
from trying to understand the recently discovered discrepancy between the critical exponents 
for matter coupled to gravity in two dimensions as computed in the lattice regularized model for 
gravity [16,17], and the corresponding conformal field theory predictions [18,19]. Particularly 
significant in this respect appears to be the recent realization that the conformal field theory 
exponents describe two-dimensional random systems in flat space, and do not correspond to 
"gravitational" dressing of correlators [21,22]. 

The plan of the paper is as follows. We shall first introduce our notation and describe 
the gravitational action, including matter fields. In order to keep our discussion as simple 
as possible, we shall limit it mostly to the two-dimensional case, although it can be easily 
generalized to higher dimensions. The subsequent sections will then be devoted to the systematic 
development of the lattice weak field expansion. The results presented here will help elucidate the 
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ABSTRACT 
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CO We develop the general formalism for performing perturbative diagrammatic expan- 

O 

\Q sions in the lattice theory of quantum gravity. The results help establish a precise 

o\ 

correspondence between continuum and lattice quantities, and should be a useful guide 
for non-perturbative studies of gravity. The Feynman rules for Regge's simplicial lattice 
formulation of gravity are then discussed in detail in two dimensions. As an application, 



the two-dimensional conformal anomaly due to a D-component scalar field is explicitly 
computed in perturbation theory. 
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